State preparation and dynamics of ultracold atoms in higher lattice orbitals 
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We report on the realization of a multi-orbital system with ultracold atoms in the excited bands 
of a 3D optical lattice by selectively controlling the band population along a given lattice direction. 
The lifetime of the atoms in the excited band is found to be considerably longer (10-100 times) 
than the characteristic time scale for inter-site tunneling, thus opening the path for orbital selective 
many-body physics with ultracold atoms. Upon exciting the atoms from an initial lowest band Mott 
insulating state to higher lying bands, we observe the dynamical emergence of coherence in ID (and 
2D), compatible with Bose-Einstein condensation to a non-zero momentum state. 
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Ultracold neutral atoms in optical lattices form a re- 
markable model system for investigations of strongly cor- 
related quantum phases, highly relevant to modern solid 
state physics [HE]. So far, experiments have mostly con- 
centrated on strongly interacting atoms in the lowest en- 
ergy band of an optical lattice. In solid state physics, 
however, orbital selective phenomena are known to play 
a crucial role for the description of topical materials 
with strong electronic correlations, such as the transition- 
metal oxides exhibiting metal-insulator transitions, su- 
perconductivity and colossal magnetoresistance. The ex- 
tension to excited Bloch bands with ultracold atoms in 
optical lattices allows to enter a regime beyond s-wave 
isotropy of the on-site atomic wave function on the one 
hand, and promises a novel route for the realization of 
long range interactions on the other hand. Recently, 
these aspects have attracted much attention and have 
lead to a series of predictions for novel effects and quan- 
tum phases that could be realized with ultracold atoms 
in higher-lying orbitals. These include the generation 
of novel multi-flavor and multi-orbital systems [3j HJ IS], 
supersolid quantum phases in cubic lattices [6j[7], quan- 
tum stripe ordering in triangular lattices |8j or Wigner 
crystallization in honeycomb lattices |9j. An essential ex- 
perimental question is however, how such systems could 
be prepared and whether they are stable. 

In this Letter, we present the experimental realization 
of a p-wave orbital system in the first excited Bloch band 
of a 3D optical lattice. By driving stimulated Raman 
transitions [10] between the different energy bands, start- 
ing from an initially Mott insulating state, we are able to 
control the population in different Bloch bands. Previ- 
ously it had been shown that for an array of quasi 2D 
degenerate Bose-gases, such excitations tend to decay 
very rapidly [llj. However, in the 3D lattice configu- 
ration, we find that such decay processes can be signifi- 
cantly suppressed, due to a reduced final density of states 
for inelastic scattering events |3j. Here, the lifetimes of 
the atoms in p-wave orbitals can extend up to hundred 
tunneling times, thus opening the avenue for orbital se- 
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FIG. 1: (a) Schematic of anisotropic tunneling following an 
excitation along the x— direction into a p-wave orbital using 
stimulated Raman transitions, (b) Rabi oscillations between 
the s— and p— wave orbital demonstrate the coherent coupling 
between the levels. 



lective physics with ultracold atoms. Starting from an 
initially Mott insulating state in the lowest energy band, 
we transfer the population in one or two dimensions into 
excited p-wave orbitals. We find that upon quenching the 
system into the superfluid regime of the excited bands, 
it exhibits a dynamical evolution to a state with long 
range coherence at non-zero quasi-momentum, providing 
a possible route for Bose-Einstein condensation to non- 
zero momentum [12J. 

We prepare our sample by first loading a quasi-pure 
BEC of typically 2 x 10 5 87 Rb atoms into a deep, sim- 
ple cubic 3D optical lattice formed by three far de- 
tuned standing wave light fields along the x—,y—, and 
z— directions. Here, the lattice potential can be approxi- 
mated by an array of separated 3D harmonic oscillators 
with discrete vibrational levels. We denote the transi- 
tion frequencies along one axis, e.g. along the x-axis, 
between the vibrational levels \n x n y n z ) and \n' x n y n z ) 
as , where n 7 is the vibrational quantum number 

along the j-axis. Population can be coherently trans- 
ferred between these levels using a stimulated two-photon 
Raman process between the initial and final external 
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states with an effective Rabi frequency = ^ A 2 • 
\{n' x n y n z \ e *( k ^i - k ^ 2 )' x | n x n y n z ) | 2 . Here, (0 2 ) de- 
notes the resonant Rabi frequency for the single-photon 
transition on the 87 Rb D 2-line to which the two counter- 
propagating Raman beams with wave vectors k^ and 
k/, 2 along the x-axis are off-resonant by the detuning A. 

For the experiments presented here, we want to restrict 
the Raman coupling to the two lowest Bloch bands along 
the x-direction only. Therefore, we initially configure the 
lattice along the orthogonal y— and z— lattice directions 
at a large lattice depth of V y = V z = bbE r , where E r 
is the single-photon recoil energy E r = h 2 /2m\f at with 
^lat = 843 nm the lattice laser wavelength and m the 
mass of a single atom. The x— axis is initially prepared 
at V x — 40 £? r , corresponding to an excitation frequency 

(x) 

of ~ 27r x 40 kHz. While loading the atoms into this 
lattice configuration, we avoid any population of higher 
bands both from heating processes during the ramp and 
from interaction-induced coupling [T3| H4] by ensuring 
that the inverse ramp time and the on-site interaction 
is much smaller than the band gap. Due to the anhar- 
monicity of the trapping potential at a lattice site, the 
excited vibrational states are not equally spaced. There- 
fore, it is possible to avoid additional coupling of the first 
to the second excited state 1 200) if the coupling strength 

(x) (x) 

Q e ff fulfills Q e ff <C c^oi — ^12 ~ 27r x 4.6 kHz for our lat- 
tice parameters. The coupling strength should however 
also be large compared to a broadening of the level spac- 
ing caused by the tunnel coupling and spatial inhomo- 
geneities due to the Gaussian beam profile of the lattice 
lasers. Then the stimulated Raman process can drive ef- 
ficient coherent Rabi-oscillations between the two states 
1 000) and 1 100) , without significant population of any 
other excited vibrational state (see Fig. [TJd) . For a beam 
waist of approximately 150 /im of the Raman laser and 
a detuning of A = 2tt x 6.8 GHz from the 87 Rb D 2 -line, 
we find typical optimal coupling strengths of ^ e ff ~ 27rx 
2.2 kHz. The Rabi oscillations are damped on a time 
scale of 2.5 ms which coincides with the calculated inho- 
mogeneity of 400 Hz due to the harmonic confinement 
across the atom cloud. From these Rabi oscillations, 
we deduce the pulse length of a Raman 7r-pulse for the 
1 000) <-> 1 100) transition for which we achieve a transfer 
efficiency of nearly 80%. 

The atomic sample is detected using absorption imag- 
ing after a ballistic expansion. When the trap is abruptly 
switched off before this time-of-flight (TOF) period, the 
momentum distribution after TOF approximatively rep- 
resents the Fourier transform of the original macroscopic 
wave function in the lattice. Alternatively, the popu- 
lation of each vibrational (orbital) band can be mea- 
sured by mapping it to the corresponding Brillouin zone 
[UJ [15j [16] when the lattice potential is ramped down 
adiabatically with respect to the on-site trapping fre- 
quency (images in Fig. [TJd). 
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FIG. 2: Measured lifetime of the excited state 1 100) with re- 
spect to (a) the total atom number for V x = 40 E r and (b) 
the lattice depth V x along the excitation axis (for 4.5 xlO 4 
atoms). The lifetime r was determined from an exponential 
decay fit to the measured population dynamics oc e _t ^ T (see 
inset) neglecting short time effects. 



The lifetime of the atoms in excited bands is a crucial 
point for the realization of many proposed novel quan- 
tum phases. In general, the population in the first ex- 
cited band constitutes a meta stable state that can decay 
by collisional events [16]. In a 3D optical lattice, how- 
ever, energy conservation combined with a small number 
of accessible final scattering states suggests that decay 
processes are suppressed. In particular, the anharmonic 
spacing of the vibrational levels should inhibit the first 
order decay process, which transfers two colliding p-wave 
bosons into one higher excited second band orbital and 
one ground state orbital. A more detailed analysis shows 
that higher order processes, i.e. involving more particles 
and higher vibrational states, can lead to a decay of the 
excited population |3j. 

In order to experimentally determine the lifetime of 
the first excited band, we measure the population in the 
different vibrational states versus the hold time in the lat- 
tice for different atom numbers and various lattice depths 
along the excitation axis. For the measurement versus 
atom number, we first prepare the system in the first ex- 
cited state 1 100) by a Raman 7r-pulse. We subsequently 
record the populations of ground and excited vibrational 
states over time. For large atom numbers, we observe an 
immediate start of the decay of population from the first 
excited band, which is well described by an exponential 
decay of the form exp(— t/r) (see inset in Fig. [2^i). For 
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decreasing atom number, the decay is delayed for short 
times. This can be understood through the atom number 
distribution in the lattice. For large atom numbers, many 
lattice sites are expected to be occupied by more than 
one atom. Consequently, the atoms can interact without 
delay. For decreasing atom number, less and less sites 
are occupied by more than one atom, until only isolated 
atoms are trapped in the lattice. In this case, atoms have 
to tunnel before they collide with another atom and sub- 
sequently decay. Excluding the short time behavior, we 
also fit the decay for small atom numbers with an expo- 
nential decay. In Fig. [2^t we plot the decay time r versus 
the initial atom number. The data shows a strong atom 
number dependence, with decreased lifetimes for larger 
atom numbers. Moreover, the data reveals a kink in the 
decay time for atom numbers of approximately 4.5 xlO 4 , 
which roughly matches the point where in our system we 
expect a large fraction of doubly occupied sites to oc- 
cur [17]. Below this atom number we observe a strongly 
rising lifetime of the system up to a factor of five com- 
pared to the higher density case, compatible with the 
interpretation that several tunneling times are required 
for collisions to occur. 

This interpretation is supported by an additional mea- 
surement of the lifetime versus lattice depth. Lower- 
ing the lattice depth implies increased tunneling rates 
which enhances the rate of atom-atom collisions, thus 
leading to reduced lifetimes. Furthermore, due to the 
larger bandwidth, the suppression of inelastic collisions, 
resulting from the anharmonicity of the potential, is re- 
duced. Fig. [2)3 shows the measured lifetime of 1 100) for 
different lattice depths along the excitation axis, normal- 
ized to the corresponding tunneling time scale h/ J' . For 
all these measurements, the sample was prepared in the 
same way at an x-axis lattice depth of 4QE r before this 
axis was ramped to the final value in 50 /is. Interestingly, 
the lifetime of the first excited Bloch band can exceed the 
time scale associated with inter-site hopping dynamics by 
a factor of almost 100. The lattice depth-dependence of 
this scaling shows a similar behavior as recently calcu- 
lated lifetimes in the first excited band, which accounts 
for for higher order decay processes [3]. 

The long lifetimes of the excited state 1 100) compared 
to the tunneling time scale should allow for the atoms to 
delocalize and exhibit spatial coherence within the first 
excited Bloch band when the system is quenched into 
the superfluid regime for this band. In order to observe 
the emergence of coherence, we excite our initially pre- 
pared Mott-insulator by a Raman 7r-pulse into the first 
excited band along the x-axis. Subsequently, we reduce 
the lattice depth in the excitation axis within 50 /is to 
a lattice depth Vhoid l ess than 25 E r . There we expect 
the lowest energy state in the first excited band to be 
superfluid. The corresponding in-trap momentum distri- 
bution is then recorded after a TOF period for varying 
hold times by switching off all trapping potentials. Right 
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FIG. 3: Emergence of ID coherence in the first excited 
Bloch band: Directly after lowering the latttice depth to 
Vhoid — 17 E r , the first excited band is still homogeneously 
populated (a) and the TOF image shows the Wannier func- 
tion of the first excited state as envelope (b). Within 1ms 
hold time, the system relaxes to a non-zero momentum state 
(d) and shows interference structure (e). (c,f) Correspond- 
ing horizontal integrated profiles with fits to the coherent con- 
tribution of the first excited state (red) and incoherent contri- 
butions of the vibrational ground state (blue) and first excited 
band (green), (g) Time evolution of the coherent and inco- 
herent contributions at Vhoid — 20 E r (color code as above). 



after lowering the lattice depth the atoms are still local- 
ized at individual lattice sites, corresponding to a homo- 
geneously filled excited band (see Fig. [3^1). In this situ- 
ation, the system features no well-defined phase relation 
between adjacent lattice sites and consequently does not 
show any interference pattern but the corresponding on- 
site momentum Wannier function of the state 1 100) (see 
Fig. j3|)) . For hold times longer than the tunneling time 
at the corresponding lattice depth, tunneling processes 
and collisions lead to the relaxation of the population 
in the first Bloch band toward the points of lowest en- 
ergy, i.e. the edges of the dispersion relation, see Fig.[S]l. 
These points correspond to a well-defined phase shift of 
Acj) = 7r between the atomic wavefunction on neighbor- 
ing lattice sites leading to the observation of interference 
fringes in TOF for odd multiples of ±ftk along the x-axis 
(see Fig. [3^), with k being the wavevector of the lattice 
laser light. Interestingly, this situation implies the relax- 
ation to a non-zero momentum state of the system fT2] . 
For hold times longer than 3 ms, the steady decay of the 
first excited vibrational level and the increasing popula- 
tion of the vibrational ground state lead to an observed 
loss of interference contrast on a time scale which corre- 
sponds to decay processes as discussed before. In order to 
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FIG. 4: Time evolution of the measured interference pattern 
emerging out of the states (a) 1 100) and (b) 1 1 10) after the x— 
and 2/-axes were ramped down to 17 Er in 50 /jls after the state 
preparation while the z-axes was kept at 55 Er. In analogy 
to the ID case, for zero hold time (a,b) the corresponding 
Wannier product wave functions of 1 100) and 1 1 10) (c+d) 
are visible, indicating an initial homogeneous occupation of 
the first excited Bloch band along (a) one and (b) two axes. 
A totally relaxed system in the two cases should yield grid- 
(e) and point-like (f) interference patterns. 



quantify our observation, we first integrate the measured 
absorption images (Figj3|),e) along the y-axis. We fit the 
resulting profiles (Fig. |3p,f ) by a model, which accounts 
for three contributions to the profile: a coherent and an 
incoherent contribution of the first excited state as well 
as an incoherent contribution of the vibrational ground 
state (see Fig. [3j:,f). For the incoherent contributions, 
we assume the corresponding Wannier functions of the 
ground and first excited states, respectively. We use a 
minimum coherence length of 6 lattice sites as a criterion 
for separating the coherent and the incoherent fractions 
of the excited band population. 

Fig. [3J5 shows the time evolution of the three contri- 
butions extracted from a fit for Vhoid — 20 E r . Within 
the first millisecond, we observe damped oscillations in 
the coherent and incoherent fractions in the excited band 
with a frequency of uo osc = 2tt x 2.6(1) kHz. The origin of 
these oscillations, which only emerge for a small range of 
the potential depths along the excitation axis, is not fully 
understood so far. Possible explanations could be inter- 
and intra-band interference effects, or the onset of an 
oscillating superfluid order parameter as predicted for a 
rapid projection of a Mott-insulating system into the su- 
perfluid regime [18j[T9]. For hold times longer than 1 ms, 
the oscillations of the coherent and incoherent contribu- 
tions are damped and reach a steady state from which 
both populations subsequently decay as discussed in the 
previous section. 

Vibrational excitations in energetically degenerate lat- 
tice axes are coupled, if the axes are not aligned perfectly 
perpendicular. This is the case in our setup. Conse- 
quently, for the same sequence as used to observe the 
emergence of coherence in ID, an excitation along 1 100) 
can couple to |010) if the potential along the y— axis is 
ramped after the Raman 7r-pulse such that V y = V x . The 



ensuing dynamics in 2D can be observed (see Fig.|4^) and 
shows an intriguing time evolution. By a slightly more 
intricate experimental sequence, the system can also be 
initially excited to 1 110) . Here, we also observe a pro- 
nounced dynamics (see Fig. |4j)) , where the momentum 
distribution in both cases finally relaxes to a grid-like 
structure after approximately 400 /is hold time. This al- 
ludes to a situation where coherence is built up along the 
two spatial directions independently. The particles thus 
seem to tunnel preferentially in one direction, rather than 
along two axes, thus yielding no emerging cross-coherence 
between the two directions. 

In conclusion, we have demonstrated the coherent ma- 
nipulation of vibrational bands in an optical lattice, 
thereby realizing a p-wave orbital system within a 3D 
optical lattice. We have measured the stability of atoms 
in the first excited Bloch band and could observe life- 
times one or two orders of magnitude larger than the 
characteristic time scales for tunneling. Finally, we have 
observed the dynamical emergence of coherence in ID as 
well as in 2D, where the system relaxes to a non-zero 
momentum state including an oscillatory dynamical evo- 
lution. Due to the long lifetimes, we believe that this 
system is a promising starting point for the realization of 
multi-orbital quantum phases with ultracold atoms. 
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